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Exact diagonalization calculations of finite antiferromagnetic spin-1 Heisenberg chains at finite 
temperatures are presented and compared to a recent inelastic neutron scattering experiment for 
temperatures T up to 7.5 times the intrachain exchange constant J. The calculations show that 
the excitations at the antiferromagnetic point q = 1 and at q = 0.5 remain resonant up to at 
least T = 2 J, confirming the recent experimental observation of resonant high-temperature domain 
wall excitations. The predicted first and second moments are in good agreement with experiment, 
except at temperatures where three-dimensional spin correlations are most important. The ratio 
of the structure factors at q = 1 and at q = 0.5 is well predicted for the paramagnetic infinite- 
temperature limit. For T < 2J, however, we found that the experimentally observed intensity is 
considerably less than predicted. This suggests that domain wall excitations on different chains 
interact up to temperatures of the order of the spin band width. 



PACS numbers: 75.10.Jm, 75.40. Gb, 75.40. Mg 
I. INTRODUCTION 

The properties of antiferromagnetic spin-1 Heisenberg 
chains result from the presence of strong quantum fluc- 
tuations in the ground state, preventing magnetic or- 
der even at zero temperature where thermal fluctua- 
tions are absent. The ground state is a macroscopic 
quantum state, a non-magnetic spin singlet, where the 
spin correlation function falls off exponentially with in- 
creasing distance.™ At zero temperature an antiferromaetj 
netic spin-1 chain has a hidden long-ranged spin ordeiou 
which can carry long-lived spin-1 triplet excitations de- 
spite the absence of conventional magnetic long-range ok-i 
der. These excitations can be regarded as domain wallsu 
of the hidden spin order and they have a mass A = 0.41J 
(J is the exchange constant) which is generated by strong 
quantum fluctuations and is not a result of an anisotropy 
in the systemE'Effi The excitation spectrum of an anti- 
ferromagnetic spin-1 Heisenberg chain is fundamentally 
different from the equivalent spin- 1/2, chain whose exci- 
tations are gapless spin-1/2 particles.Q 

At finite temperature, temperature fluctuations dis- 
order the hidden "quantum order" and this leads to a 
renormalization of the triplet excitations. For T < A, 
the triplet excitations are well-defined and the spin- wave 
velocity is constant, but their mass - defined as the exck 
tation energy at the antiferromagnetic point - increases .u 
This upward renormalization of the excitation energy is 
well describedJpy_a low-energy field theory, the non-linear 
sigma model.Qa'EHl For finite temperatures, the excita- 
tions acquire a finite life-time due to collisions with ther- 
mally excited particles and the temperature dependence 
of the decreasing life-time of the excitations is well de- 
scribed by a semi-classical model .of .&ai>ped spin-1 parti- 
cles with a relativistic dispersionOt§E3 

For T > A, there is a cross-over to a qualitatively 
different behavior. In this temperature range.,the exci- 
tations are resonant at least up to T = 2.7J,0 but the 



excitation width is nearly as large as the excitation en- 
ergy and the spin- wave velocity decreases with increasing 
temperature. These high-temperature excitations can-|be. 
regarded as domain walls of a local hidden spin ordenl3'E2l 
which is the remainder of the hidden long-range order at 
T = K. However, a quantitative theoretical under- 
standing of the high-temperature spin dynamics is miss- 
ing at present. 

In the lack of reliable analytical methods for a quanti- 
tative calculation of dynamical properties at finite T > A 
we calculated the scattering function numerically by ex- 
act diagonalization (ED). We considered finite spin-1 
chains with an antiferromagnetic Heisenberg exchange, 
described by the Hamiltonian 

H = J^2Si-Si +1 , (I) 

i 

and calculated the dynamic structure factor and fre- 
quency moments for temperatures up to T = 7.5 J. The 
results of the calculations were then compared with the 
measured neutron scattering spectrum of the antiferro- 
magnetic spin-1 Heisenberg chain compound CsNiCbj. 

CsNiCbj crystallizes in a hexagonal crystal structure, 
Dgt space group and is one of the best model systems for 
investigating antiferromagnetic spin-1 Heisenberg chains. 
The spin-1 components are carried by Ni +2 which inter- 
act via a super-exchange involving Cl _ -ions. The super- 
exchange interaction J along the c-axis is much stronger 
than the super-exchange interaction J' in the basal plane. 
Thus CSN1CI3 approximates a system of weakly coupled 
spin-1 chains. The Hamiltonian can be written as 

i 

+ J'£S i -S J -£>£ 4 (S?) 2 , (2) 

<i,j> 

with J = 2.28 mcV and J' = 0.044 mevH&El The 
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sum in the basal plane includes the interaction between 
nearest- neighbors < i,j > only. The single- ion crystal- 
field anisotropy D — 4 /ieV is small enough that the 
exchange interaction between the spins is isotropic in spin 
space. A detailed description pi-Uac. neutron experiments 
has been presented elsewhere. ErE3ll3'E£l 

This paper is structured as follows: Section II intro- 
duces the numerical calculations, and shows that the ex- 
citations are resonant up to at least T = 2 J. In Section 
III energy-integrated quantities of the dynamic structure 
factor are quantitatively compared with the experiment 
and effects of biquadratic and next-nearest neighbor in- 
teractions on the excitations are discussed. Section IV 
addresses the temperature dependence of the expectation 
value of the Hamiltonian and shows that there is excel- 
lent agreement between experiment and theory. Section 
V summarizes the results. 



II. EXACT DIAGONALIZATION 

The scattering function was obtained numerically by 
ED for single spin-1 chains of finite length through the 
following formulae: 
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where \i) are the exact eigenstates of the Hamiltonian, 
with energy u>i, and Boltzmann population factors rii. 
N s is the number of spins in the chain, and S z (R) is the 
operator for the z component of the spin at position R. 
q is the wave-vector along the chain. 

The calculation of finite temperature properties re- 
quires knowledge of the full spectrum of the Hamilto- 
nian so that only rather small sizes can be handled. We 
used a chain with N s = 8 sites with periodic bound- 
ary conditions. The dimension of the Hilbert space is 
3 s = 6561. To reduce the numerical effort, we exploit 
the conservation of the z component of the total spin 
gtot _ S Z (R), which makes the Hamiltonian matrix 
block-diagonal in the basis of eigenstates of Sl ot . Since 
the operator S z (q) commutes with S^ ot , matrix elements 
in Eq. ^ are nonzero only if \i > and \j > have the same 
value of Sl at . Thus, once the partition function (enter- 
ing rii and rij) has been determined, the contribution of 
each S"* * eigenspace to the scattering function can be 
calculated independently. 
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FIG. 1. The solid line represents the calculated dynamic 
structure factor S zz (q,ui) for q = 0.5 (upper figure) and 1 
(lower figure) at T — 2 J as a function of energy u)/J (top 
axes). The calculations for a finite chain of 8 sites were con- 
voluted with a Lorentzian with a half width 0.3J in order to 
smooth out finite-size oscillations. The intensity is for a chain 
ofiV s =8 spins and such that J2k- J dcjS zz {ki,iv) = s(s + 1)JVs 
for S = 1 (right axes). The measured neutron scattering spec- 
tra are shown as solid circles and as a function of neutron 
energy transfer uo (bottom axes). The data were measured 
at 50 K ~ 2J usuae-ithe RITA and the MARI spectrometer, 
taken from Refs.B'E3, and in units shown on the left axes. 
The axes on the left and right in the two figures were offset 
to account for the experimental background scattering. The 
dashed lines corresponds to the elastic scattering peak (not 
completely shown) which arises due to incoherent scattering 
at the sample and its environment. 



Zero-temperature results have been obtained for sizes 
up to N s — 18 by using the Lanczos algorithm to find 
the exact ground state of the system. In these calcu- 
lations, the Hamiltonian matrix is further reduced to 
block-diagonal form by exploiting translational invari- 
ance, using a basis of states with given total wave- 
vector K tot = h = 0, 2n/N s , ....,2n(N s -l)/N s . Then, 
(i\S z (q)\j) ^ only if S*°»(t) = S\°\j) and K tot {i) = 
K~tot(j) + q- The T = dynamical structure factor can 
be obtained without performing an explicit calculation 
of all the spectrum of the Hamiltonian, by using in place 
of Eq. (1) a continucd-iraction representation of the dy- 
namical susceptibilityLH, whose coefficients can be calcu- 
lated recursively once the ground state is known. We 
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verified as a check of the caiculations that the sum rule 
jdu)S zz (q,uj) = (S z (-q)S z (q)) is verified. 

Fig. [j] shows the calculated S zz (q,u>) as a function of 
energy transfer ui/J for q = 0.5 and q = 1 and for a 
temperature T = 2 J. The energy spectra were obtained 
by convoluting the discrete spectra for finite chains with 
Lorentzian functions (width 0.3 J), q is given in units 
such that q — 1 corresponds to the antiferromagnetic 
point, which in theoretical work is often referred to as 
the 7r-point. Throughout the paper, the intensity of the 
calculated structure factor S zz is given for a chain of 
N s spins and normalized so that Y]u. J dujS zz (ki,uj) — 
s(s+i)N s because m an isotropic system one spin com- 
ponent carries exactly one third of the intensity. 

The calculated excitation spectra are broad but still 
resonant at T — 2J, consistent with neutron scattering 
measurements which showed that the Haldane excitation 
survives as a resonant feature up to at least T = 2.7 jBt£l 
The calculated scattering weight at q = 1 is centered at 
1.5 J so that the [Scattering is at higher energies than at 
low temperature.EI At q = 0.5, the increase of the temper- 
ature has the opposite effect and the calculated scattering 
is centered clearly below the low-temperature maximum 
of the dispersion which is uj — 2.71J. 

The calculated excitation spectra are compared in 
Fig. [j] to neutron scattering data from CsNiCl3.Erll3 The 
chains in CsNiCi3 run along the c-axis, so the third com- 
ponent I of the wave-vector transfer Q = (h, k, I) corre- 
sponds to the wave- vector transfer along the chain. For 
the remainder of this paper q denotes both the wave- 
vector transfer along a single chain used in the theoret- 
ical calculations and the wave-vector transfer along the 
chain in CsNiCbj, thus I = q. The data in Fig. [j] was 
measured at T = 50 K which corresponds to T = 2 J. 

Three-dimensional (3D) spin-spin interactions - due 
to the interchain coupling J' - are most important for 
temperatures T < 15 K and they induce aatiferromag- 
netic long-range order below Tn = 4.84 K.L3 Above Tn, 
3D interactions lead to a dispersion.-Qf the excitations 
perpendicular to the chain directionJiJ which decreases 
with increasing temperature but 3D spin correlations 
may be still be important at T — 50 K. The measure- 
ments for q = 1 were performed at a wave-vector trans- 
fer Q 1D = (0.81,0.81,1) where th&,Fourier transform of 
the interchain coupling vanishesuBj and where, within a 
Random Phase Approximation, (RPA), the chains behave 
as if they were not coupled.Erfl So to first order the ex- 
periments probed the dynamics of isolated chains. For 
q = 0.5, the interchain coupling has only a negligible 
effect on the dynamic structure factor. Thus it is not im- 
portant that the wave-vector transfer perpendicular to 
the chain varies in energy scans with the MARI spec- 
trometer and the observed scattering for q = 0.5 can be 
directly compared to the calculations. 

As shown in Fig. [j] the calculated and observed dy- 
namic structure factor are in good agreement for q = 0.5 
and q = 1 and the calculated excitation energies and 



widths are well predicted by theory. The small differences 
between experiment and theory may arise from finite-size 
effects in the calculations. 
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FIG. 2. Ro(q) = S y >0 (g), Ri(q) =< lo >„>„ (q) and 
R2(q) = \/ < u' 2 >ui>o(q) for q — 1 obtained from ED calcu- 
lations at zero temperature as a function of the inverse chain 
length 1/N S . The solid lines are fits to the data to extrapolate 
to an infinite chain as explained in the text. 

The positive frequency moments of the spectrum are 
expected to be more reliable than the line-shape of 
S zz (q,u). For example, at T = the ratio of the con- 
tinuum to the Haldane mode contributions for a single 
chain is almost independent on the chain length and is 
well reproduced in small chain calculations, despite the 
fact that the line-shape of the continuum in S zz (q = 1, u>) 
is not correct. The following moments were considered: 

Ml) = S zz >0 (q) = J u>0 S zz {q,Uj)d^ , (6) 

Rt(q) =< uj > ul>0 (q) = 
Sg >0 (q)- 1 f u>0 wS"(q,u)dLJ, (7) 



Ra(q) = V< >^>o(q) = 

(stub)' 1 f»>o^ st 'b> u ) du ) 1/2 • ( 8 ) 

For k-eT > 0.5J finite size effects in Ri appear to be 
small (this is seen by comparing results for N s = 4, 6, 
8). For smaller values of T, finite-size effects appear to 
be sizeable for q = 1 and very small for q = 0.5. We 
extracted the values of Ri(q = 1) at T = K for N s = oo 
by performing a second-order finite-size scaling fit: 

R\(N S ) = Ri (N s = oo) + Ai/N s + BjNl (9) 

Ai and Bi were calculated from Ri for N s = 14, 16, 18. 
The extrapolated Ri(N s ) were compared with the cal- 
culated Ri(N s ) for N s — 6, 8, 10, 12 and the results in 
Fig. H show that Eq. [)] holds down to rather small sizes, 
with just a small deviation for Ro(q)(N s = 6) (associated 
with a small 1 /Nf correction) . 
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III. QUANTITATIVE COMPARISON WITH 
EXPERIMENT 

A. Excitations at q = 0.5 

For q — 0.5, the coupling of the chains in CSMCI3 has 
a negligible effect on the dynamic structure factor and a 
direct comparison between the experiment and the ED 
calculations is possible. Fig. [3] shows the experimentally 
observed Ro{q), Ri{q) and i? 2 (<z) for q = 0.5 for tem- 
peratures between T = 6 and 200 K. The data were de- 
termined numerically from the neutron spectra measured 
using the MARI spectrometer after subtracting the flat 
background and accounting for theiaiiagnetic form factor 
of the neutron scattering intensity.ll2l 
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FIG. 3. (a) Integrated neutron scattering intensity Ro(q) 
for positive energy transfers at q — 0.5 as a function of tem- 
perature on a semi-logarithmic plot. The data were mea- 
sured using MARI and for [temperatures between T = 6 and 
200 K and taken from Refill The intensity was determined 
by adding up the observed spectra for positive energy trans- 
fers after subtracting theJlat background and accounting for 
the magnetic form factorJij The solid line corresponds to the 
predicted Ro(q = 0.5) obtained from ED calculations and fit- 
ted to the data, and it is shown in units shown of the right 
axis, (b) First energy moment Ri(q = 0.5) of the neutron 
scattering spectra observed on the energy loss side, (c) Sec- 
ond energy moment J?2(g = 0.5) of the scattering spectra on 
the energy loss side. Solid lines in (b-c) are ED predictions. 



The ED results for R~o(q — 0.5) were adjusted in a 
least-squares fit to the experimentally observed intensi- 
ties using an overall scaling factor. No scaling factors 
were used for R\{q — 0.5) and R 2 {q = 0.5) which scale 
with the intrachain exchange J = 2.28 meV, accurately 
known from high-field magnetization measurements. 

The agreement between theory and experiment is ex- 
cellent as shown in Fig. |[ The ED calculations repro- 
duce the flat temperature dependence of the intensity 
at q = 0.5. A comparison with the overall scale is not 
possible because S(Q,uj) was not measured in absolute 
units. Excellent quantitative agreement is found for the 
temperature dependence of the first and second moment 
at q = 0.5, as shown in Fig. ||(b)-(c). The small dif- 
ference between the experiment and the calculations for 
T < 10 K, where R x {q = 0.5) and R 2 (q = 0.5) are slightly 
lower than predicted by our ED calculations, is at least 
partly due to finite-size effects: our ED calculations of fi- 
nite chains overestimate the energy oj(q = 0.5) = 2.77 J = 
6.31 meV, whereas a more accurate zero-tempeca.ture cal- 
culations give uj{q = 0.5) = 2.71J = 6.17 meV.ll3 

B. Excitations at q = 1 

Fig. |] shows Ro(q), Ri(q) and R 2 (q) for q = 1, mea- 
sured between T — 6 and 70 K using the RITA spectrom- 
eter and compared to ED calculations of a finite chain of 
8 sites. The neutron scattering measurements.-iYere per- 
formed at the non-interacting wave- vector Q 1D &E£I where 
the chains behave within RPA as if they were not coupled. 
To first order, the temperature dependence of Ro(q = 1), 
Ri(q = 1) and R 2 (q = 1) are thus expected to match 
that of a single chain. 

Ri(q = 1) is about 0.65 J at zero temperature, and it 
is almost three times as large at T = 2J. This reflects 
the rapid increase of the, excitation energy at q = 1 with 
increasing temperatureJa The zero-temperature value of 
Ri(q = 1) is higher than the Haldane gap 0.41J due to 
the presence of a weak sinderchain multi-particle con- 
tinuum at higher energies There is good agreement 
between the measured and calculated Ri{q = 1) and 
R 2 (q = 1) for temperatures T > 15 K as shown in 
Fig. ^. Below this temperature, the experimentally ob- 
served i?i(<z = 1) and R 2 (q = 1) are consistently higher 
then predicted by the calculations. In this temperature 
regime, finite-size effects are important which make the 
calculated Ri(q — 1) and R 2 (q — 1) (see Fig. ||) larger 
than those of an infinite chain and finite-size effects can 
thus not explain the difference between the experiment 
and the calculations. 

The difference between the experiment and the ED cal- 
culations are probably due to interchain couplings, which 
are known to have a strong effect on the intensity of 
the excitations just above the ordering temperature, as 
shown in Fig. 5 of our experimental paper.Ej Also, for 
temperatures T < 12 K a multi-particle continuum was 
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observed, which increases Ri(q = 1) and i?2(<Z = 1) 
due to the presence of scattering weight at high energy 
transfers. 




FIG. 4. (a) Integrated neutron scattering intensity 
Ro(q = 1) observed at the non-interacting wave- vector 
Qid = (Oj&i, 0.81, 1) using the RITA spectrometer and taken 
from Ref.£3. The solid line corresponds to the calculated 
Ro(q = 1) obtained by ED calculations. The prediction was 
scaled to match the intensity at T = 70 K and is shown in 
units shown on the right axis. The error bars are smaller 
than the symbol size, (b-c) First and second energy moment 
Ri(q — 1) and R2(q ~ 1) for the energy loss neutron scatter- 
ing spectrum observed at Q 1D using the RITA spectrometer. 
The solid line was obtained by ED calculations and the tem- 
perature and energy scaling was done using the intra-chain 
exchange J = 2.28 meV known from high-field magnetization 
measurements. 



Figs. |^(a) and || show the experimentally observed 
Ra(q) at q = 1 as a function of temperature. The 
data were measured using the RITA triple-axiSj-and the 
MARI time-of-flight spectrometer, respectively.BEJ The 
ED calculations predict that at T = 7. 5 J the ratio 
Ro(q = l)/Ro(q = 0.5) is equal to 1.26 and this can di- 
rectly be tested with the neutron scattering experiments. 
At T = 7.5J ~ 200 K the MARI measurements show 
that Ro(q = l)/Ro(q = 0.5) = 1.5(4), in good agreement 
with the calculations. At lower temperatures, however, 



the experimentally observed Ro(q = l)/Ro(q = 0.5) is 
less than predicted and this is shown in Fig. [5| Since the 
temperature dependence of Ro(q = 0.5) is in agreement 
with the calculations (Fig. ||), we conclude that the dis- 
agreement results from a reduced intensity Ro(q) at the 
antiferromagnetic point q = 1. The disagreement is most 
pronounced for T < 50 K and may be due to the fact 
that the measurements were performed close to but not 
exactly at the non-interacting wave- vector Qid- 
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FIG. 5. Neutron scattering intensity at q = 1 observed us- 
ing MARI for temperatures between T = 6 and 200 K.as solid 
circles on a semi-logarithmic plot and taken from Ref.tJ. The 
intensity was determined by integrating the observed spec- 
tra for positive energy transfers after subtracting a flat-back- 
ground and correcting for the magnetic form factor£j The 
solid line is the ED result for Ro(q) at q — 1. The scaling was 
chosen to match the prediction Ro(q = l)/Ro(q = 0.5) and is 
thus in the same units as the scaling in Fig. |3[ 



Similar result were obtained, however, from an analy- 
sis of the RITA measurements which were performed ex- 
actly at the non-interacting ID point Qid m reciprocal 
space so that to first order the neutron spectrum should 
not be influenced by 3D interactions in the system. For 
the RITA measurements the scaling factor for Ro(q = 1) 
can only be guessed because reliable measurements at 
q = 0.5 were not performed. Choosing the scaling factor 
such that the ED results match the observed Ro(q = 1) 
at T = 70 K - which is known to be approximately true 
from the MARI measurements - the solid line shown in 
Fig. ||(a) is obtained. The discrepancy between the ex- 
periment and the numerical calculations is similar to that 
observed using the MARI spectrometer. At low tempera- 
tures there is a difference between theory and experiment 
of a factor of 1.6. 

Finite-size effects may be important for T < 15 K, 
but these would underestimate the measured intensity 
rather than overestimate it (Fig. ||) and so they can- 
not explain the discrepancy between the experiment and 
the numerical calculations. This means the difference 
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between the experiment would be even bigger if longer 
chains had been considered in the calculations and the 
observed Ro(q = 1) is at least a factor of 1.6 less than 
that calculated. 



C. The influence of next-nearest neighbor, 
biquadratic and interchain spin-spin interactions 

We examined the effect of next-nearest neighbor J nn 
and biquadratic Jbi q exchange interactions on the dy- 
namic structure factor at zero temperature. The maxi- 
mum in the dispersion u>(q = 0.5) was calculated numer- 
ically for a spin-1 Hciscnberg chain with Hamiltonian 7i 
(Eq. |l|) and with additional terms 



Jnn ' S'i+2 



for next-nearest neighbor interactions or 



T~tbiq — Jbiq ^^(^i ' ^i+l)^ 



(10) 



(11) 



is reduced with respect to the value for uncoupled chains 
if the g-integrated intensity of S(q) is to be unchanged. 
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for biquadratic spin interactions. The calculations were 
performed for chains with 12 sites. The calculated en- 
ergies were then compared with the experimentally ob- 
served energies for q = 0.5 at low temperatures in order 
to estimate the strengths of next-nearest neighbor and 
biquadratic exchange interactions. 

uj{q — 0.5) decreases for increasing J nn / 'J and Jbiq /J 
and increases when these exchange couplings are nega- 
tive, as shown in Fig. ^|(a). The measured values for 
uj(q = 0.5) are between 5.9 and 6.2 me|£, foe-various ex- 
periments with 6 K < T < 12 K (Ref.BEJO). This is 
consistent with < J nn / J, Jbiq/ J < 0.05 and this range 
is shown in Fig. ^(a)-(c) as shaded area. The associ- 
ated decrease in Ra(q = 1) is at most a factor of 1.05 
(Fig. |6|(b)), which is much smaller than the observed 
reduction of Ro(q = 1) by a factor of 1.6 at low tem- 
peratures. This shows that next-nearest neighbor and 
biquadratic exchange interactions cannot lead to the ex- 
perimentally observed decrease in Ro(q = 1). 

Fig. |s|(c) shows—that the recently observed multi- 
particle continuumtJ at high energy transfer carrying 
12% of the total scattering at q = 1 cannot be caused by 
biquadratic or next-nearest neighbor interactions. This 
is because for < J nn /J Jbiq/ J < 0.05 the strength of 
the continuum is never more than 3% of the total scat- 
tering at q = 1. We note that our results for the valence 
bond solid model (Jbiq/ J = 0.33) are consistent with a 
previous study.Ea 

At low temperatures T < 12 K, the observed reduction 
of Ro{q = 1) for Q 1D is related to the correlation length £ 
of the coherent state in the chains. It was shown that in 
CsNiCU-J is reduced from the value expected for single 
chains,t3 which results in a wider peak of the structure 
factor S(q) centered at q = 1. This means that S(Q 1 - D ) 
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FIG. 6. (a) Zero-temperature ED results for the excitation 
energy at q — 0.5 as function of J nn /J and Jbiq/ J. (b) ED 
results for the integrated intensity at q = 1 and at q — 0.5 
at zero temperature as a function of Jnn/J and Jbiq/ J. (c) 
Zero-temperature results for the ratio of continuum intensity 
■fcont to total intensity Ttot at q — 1 as a function of J nn / J and 
Jbiq/ J- The shaded area in (a)-(c) correspond to the possible 
Jnn/J and Jbiq/ J range for CsNiCl3. 

The reduction of Ro(q = 1) is also related, to the 
multi-particle continuum at higher energies The 
Hohenberg-Brinkman sum rule predicts that the first en- 
ergy moment F(q) = J duj u S(q,u>) of one-dimensional 
magnetSp-with nearest-neighbor exchange interaction be- 
haves ascil 

F(q) = ~<H>(l-cos(Trq)). (12) 

It will be shown in the Appendix that this relation also 
holds for coupled chains with small interchain interac- 
tions. If the integrated intensity for a particular q is 
reduced, intensity necessarily has to be transferred to 
high energies so that the spectrum is consistent with the 
Hohenberg-Brinkman sum rule. This is because intensity 
at high energies contributes more to the first moment 
than scattering at low energies. 

The comparison between the experiment and numeri- 
cal calculations suggests that 3D correlations lower the 
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intensity at the non-interacting ID point Q 1D for T < 2 J. 
Conservation of the total magnetic intensity then implies 
that this intensity is transferred to other parts of recip- 
rocal space. The 3D dependence of the magnetic scat- 
tering suggests a coupling of domain walls on different 
chains for temperatures up to at least T — 2 J. Fur- 
ther the coupling of the chains in antiferromagnetic spin- 
1 Heisenberg chains has a larger effect on the excitation 
spectrum and extends to higher temperatures than as- 
sumed so far. Possibly the hidden string order plays a 
decisive part in this effect because it allows the existence 
of spatially-extended coherent states to higher tempera- 
tures than would be possible in magnets without a similar 
spin string order. 

IV. TEMPERATURE DEPENDENCE OF < H > 

In our experimental paper0 the first energy moment 
sum rule was used to determine the temperature depen- 
dence of the expectation value of the Hamiltonian < Ti > 
of a single spin-1 chain. Here we show that the first mo- 
ment sum rule for single chains is also valid for weakly 
coupled chains. We also compare the experimentally de- 
termined values for < Ti > to ED results. 

The first energy moment of CSN1CI3 based on the 
Hamiltonian given in Eq. [2] can be written as (see the 
Appendix) 

F(Q) = -AJF{\ - cos(Q ■ c/2)) 

- 4J'F'^(1 - cos(Q • Oj)) , (13) 
3 

where a.j are the lattice positions of the nearest neighbors 
and 

F = F a = ^2< S™.S?. +c/2 >, 

i 

i 

for a = x, y, z. The second term in Eq. |l3|is always much 
smaller than the first term, mainly because J' = 0.015 J, 
but also because Max{^ .(1 — cos(Q • a,))} = 4.5 is only 
2.25 times Max{(l - cos(Q • c/2))} = 2. In addition, 
the wave-vector transfer in the experiments was mainly 
along the c-axis, which reduces the second term. 
The dominant term in Eq. [l3] is therefore 

F(Q) = -4JF(1 - cos(Q • c/2)) 

= -4JF(1 - cos(7r g )) , (14) 

which is directly proportional to < Ti >= 3JF with Ti 
given in Eq. [j]. The expectation value of the Hamiltonian 
for a single chain can thus be obtained by fitting Eq. [l4| 
to the g-dependence aLthe experimentally observed first 
energy moment F(q)i3 
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FIG. 7. The experimentally observed < Ti > = 
J J^. < SiSi+i > as solid circles and as a function of-tempera- 
ture on a semi-logarithmic plot and taken from Ref.u. <TL> 
was determined from the first energy moment measured at the 
different temperatures. The error bars are smaller than the 
symbol size. The solid line is the ED result for the tempera- 
ture dependence of the mean energy per spin <TL> /N 3 . 

Fig. [7] ahpws the experimentally determined < Ti. > 
from Rcf.EEl, compared with our ED calculations. We 
find that there is excellent agreement between theory and 
experiment for the whole temperature range. A scaling 
factor was used to reach the agreement since < H > was 
not measured in absolute units. 

The calculated ground state energy per spin at zero 
temperature, < H > /N s = —1.41(1) J, is consistent 
with a previous numerical study of much longer chains.u 

V. CONCLUSIONS 

We have performed exact diagonalization calculations 
of a finite antiferromagnetic spin-1 chain with 8 sites. 
The calculations confirm the recent experimental obser- 
vation that the excitations are resonant for temperatures 
up to at least T = 2J. The temperature dependence 
of the calculated first and second energy moments for 
positive energy transfers is is-good agreement with the 
measured neutron scatteringpl except at low tempera- 
tures where spin correlations between chains are most 
important. Here the measured first and second moment 
for q = 1 are higher than predicted due to the presence 
of a multi-particle continuum at high energy transfers. 

The measured intensity of the excitations at the non- 
interacting wave-vector is much weaker than that calcu- 
lated for isolated chains, particularly for T < 2J. This 
suggests that solitons on different chains interact for tem- 
peratures up to at least T — 2 J. The experimentally ob- 
served < H > is in excellent agreement with exact diag- 
onalization calculations for the entire temperature range 
considered in this study (0.25 J < T < 7. 5 J). 



7 



APPENDIX 



The first energy moment is defined as 
F aa (Q) — T^ojS aa (Q,uj) — 

~ m < [[ n ' S o\ > S -q] >> 



(15) 



where the Fourier transform of the spin operators is 
defined as 



S%=Nr 1/a 52eKp(-iQri)S? t 



(16) 



Due to the discrete rotational symmetry D\ h in the mag- 
netically disordered phase, spin correlations along all 
three equivalent direction in the basal plane are the same 
and 



FL = F: 



3 = 1,2,3. 



For an isotropic spin chain J a — J, F a — F and F' a = F' 
for a = x, y, z and the first moment of S(Q, uS) is 



F(Q) = -4JF(1 - cos(Q • c/2)) 



4J'F'^(1 - cos(Q • aj)) . 



(19) 



In the the case of CsNiCl3 , the Hamiltonian includes both 
intrachain and interchain spin interaction, and it can be 
written as 

+ J x s ri s rt+a ^ + JySy.s^. i+a . + J z s r .s ri+a ^ . 

(17) 

In a cartesian coordinate system, the positions of the 
neighbors along the three directions in the basal plane 
are given as 

ai = (a, 0,0), 
a 2 = (a/2, VSo/2,0), 
a 3 = (-a/2,\/3a/2,0) . 

The first moment for S XX (Q, u>) can be split up into two 
terms, one containing the intrachain interactions J a and 
the other the inter-chain interaction J' a : 

F XX {Q) = -Uy{F y - F z cos(Q ■ c/2)) 
-2J z (F z ~F y cos(Q-c/2)) 

-Y.\ 2J y( F L 3 -Ki C0 <Q- a i)) 



+2J' z (F^-F yd co S (Q- aj ))] 



(18) 



F a and F' a • are the expectation values for spin-spin cor- 
relations along the chain direction and along either of the 
three direction in the basal plane: 

i 
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